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Abstract 



In this paper we analytically study the pricing of the arithmetically averaged Asian 
option in the path integral formalism. By a trick about the Dirac delta function, the 
measure of the path integral is defined by an effective action whose potential term is 
an exponential function, i.e. the Liouville Hamiltonian, which can be explicitly solved. 
After working out some auxiliary integrations involving Bessel and Whittaker functions, 
we arrive at the spectral expansion expression of the value of an Asian option. 



1 Introduction 



It has been known for a long time that the path integral formalism can be applied to 
the pricing of financial securities. The more standard methods in quantitative finance are 
the stochastic calculus and partial differential equations. In almost all cases the relevant 
equations are diffusion type, whose solution is determined by the heat kernel. It is well- 
known that [J] the heat kernel can be written in terms of a path integral. This is the 
starting point for the financial applications of the path integral formalism. In [2j |3] it has 
been applied to the European options and the one-factor term-structure models. In [U [5] 
it is shown how the pricing of path-dependent options can be incorporated into the path 
integral formulation. It is applied in [6] to models with stochastic volatility, and in [7] 
by the same author to the Heath- Jarrow-Morton model of forward interest rates. See e.g. 
[8]- [15] for more works in this direction. 

Among many exotic options in the financial market, the asian option is a very popular 
one. Its payoff depends on the arithmetic average of the price of the underling asset 
during the life of this option contract. The Asian options has the advantage that it is 
usually less expensive than standard options due to its smaller volatility, and its value is 
harder to be manipulated by a large market participant. So it is more safe to hold it. 
On the theoretical side, the exact pricing of the arithmetically averaged Asian option is 
a challenging problem, since the arithmetical average of a stochastic variable, which is 
logarithmic- normal, is not logarithmic-normally distributed any more. In the pioneering 
work [16J, Geman and Yor derived a closed form expression for the Laplace transformation 
of the value of the Asian option. In [17J Linetsky obtained a spectral expansion expression 
of its value in terms of confluent hypergeometric functions. This is the main motivation 
of the present paper. 

In this paper we will study the the arithmetically averaged Asian option in the path 
integral formalism. In section 2, we review the path integral formulation for a gen- 
eral path-dependent option. For the Asian option, the resulting effective action is the 
(imaginary-time) Liouville quantum mechanics. In section 3, we obtain the correspond- 
ing heat kernel by solving the differential equation it satisfies. We work out two auxiliary 
integrals in the following section 4 and 5. In section 6, we specify the payoff function and 
study the value of the put and call options. For the put option we can directly calculate 
its value, which is equivalent to the result of [17] . while for the call option we have to 
reply on the put-call parity relation. 
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2 Path integral formulation 



If S denotes the price of a stock, it is commonly assumed that X := log 5 follows the 
rule of a Brown motion. A stock option is a kind of financial derivative whose value O 
depends on the behavior of the underlining stock price [18]. For a general path-dependent 
option, its payoff function is a functional of X. That is the final yield of this option 
contract depends on the whole history of the stock price before its maturity. The usual 
European option is just a special case, whose payoff functional is local, i.e. only depends 
on the stock price at the maturity day. 

The value of a (newly written) stock option can be written in a path integral form as 



/oo pX(t)=x j i rt /jy \ 2 1 

. Ji 'L MH T^lt + rr w (21) 

In the above equation, r denoted the risk-free interest rate, a is the volatility of the stock 
price, and \i := r — er 2 /2. All of these parameters are assumed to be constant. The 
meaning of t is the remaining time of this option contract, with t = denote the maturity 
day. By introducing the following dimensionless combinations 

R=- 2 , r = *H, u=%, (2.2) 
we can simplify the above equation as 

/oo pX(t)=x f 1 PT ~\ 
dx ' e u( X '-x)/2-^r/8 / DX exp I / X2 dr \ $ ^ _ ( 2 3 ) 

-oo JX(0)=x' I 2 Jo J 

If the payoff functional takes the form $>[X] = <p (V[X]) with0 

V = V[X] = [ e x dr , (2.4) 
Jo 

then we have 

T Jo 

= - / rf(^)x- / e^ v dq. (2.5) 

T J0 Am Je-ioa 



We can restrict the integration range to (0, oo) in the first line because the functional 
V = V[X] is always positive. In the second line we have use the Laplace transformation 



J In this paper we will call this kind of options as arithmetically averaged Asian- type option or simply 
Asian option. In principle our formalism applies for any payoff function (j>. 
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of the Dirac delta function with e being any positive real number. Inserting (I2.5P into 
(12.31) we have 



— Rt roo poo re+ioo 

* KlT J-oo JO Je-ioo 

X(t)=x (■ j-T /i 



x Mexp - / l-X 2 + qe x )dr\. (2.6) 

Jx{0)=x> I Jo V 2 /J 

This is our path integral formulation of the valuation of the Asian option. We see that 
the system is driven by an effective action 

A eS [X] = £ Qx 2 + ge x ) dr , (2.7) 

which is called (imaginary-time) Liouville quantum mechanics [T9] . 



3 Heat kernel 

Define the heat kernel 



IC(T,x,x';q) := [ DXexp\-[ ( \x 2 + q eA dr\ . (3.1) 

JX{0)=x' I JO V 2 J J 

By use of Feynmann-Kac theorem, it satisfies the following initial value problem 

_dlC _ eX)C 
dr 2 dx 2 ^ 6 

fC\ T =Q = 8(x — x') . (3.2) 

Actually this is the (imaginary-time) Schrodinger equation of the effective action ( 12.71) . 
We may use the method of spectral expansion to construct the heat kernel K.(t,x,x'; q). 
Firstly we solve the following eigenvalue/eigenf unction problem 

For q G C with |argg| < it, the normalized eigenfunction is 

ifjJx) = -\/usmh(iru) K iu {\/%qe xl2 ) , u>0, (3.4) 

where K iu is the modified Bessel function of the second kind. Then the heat kernel can 
be written as 



JC(r,x,x';q) = / e~ u T / 8 if> u (x) ^ u {x') du 
Jo 

1 f°° 

= — e- u2r/8 K iu (^e x/2 )K iu (^e x ' /2 )smh(iru)udu. (3.5) 
71 Jo 
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It can be explicitly proved that, when q > 0, (13. 5p is indeed the solution of (13. 2p . That 
it satisfies the differential equation can be easily checked. In the appendix we will show 
that it also satisfies the initial condition, i.e. the completeness of {ipi U {x) \u > 0}. Due 
to the inverse Laplace transformation in (I2.6p . we need q to be complex with Reg > 0. 
Since the solution of (13. 2p should be a holomorphic function in the complex g-plane cut 
open along the negative real axis, we can insert (I3.5P into (12.61) to calculate the option 
value, and the result turns out to be correct. Therefore the expression of the option value 
becomes 

—Rt roo i>oo pt+ioo 

0{t,x) = —— e~ ux /2 ~ u2t/8 / dx'e ux ' /2 ^0(0/ d Q ^ 

%KIT J-oo Jq Je-ioo 

x— / e- u2r/8 K iu (J8q'e x/2 )K iu (y8q'e x ' /2 )smh(TTu)udu. (3.6) 
^ Jo 



4 Integrating out x' 



In this section we will consider the integration over the variable x' . Define 



M(r,x;q) :~- 



3 ra ' /2 /C(r, x, x'\ q) dx' , 



(4.1) 



which satisfies the following initial value problem 

DM 1 d 2 M 



dr 

M \ T=0 



2 dx 2 

ux/2 



+ qe x M 



(4.2) 



To solve this problem we use the following expansion of the initial configuration [2T 

2 



,ux/2 



7T 2 q v l 2 



l>CO 






r (^) 


Jo 





Ki U (\/8q e X//2 ) sinh(7ra) u du 



+ 



2l-f/2 



7 u/2 



[-I//2] 

E 

n=0 



'-v-2n) 



n\ r(— v — n + 1 



K- u - 2n (V8qe 



x/2\ 



(4.3) 



Note that when v is not positive, the function e ux l 2 is not in L 2 (IR), so the above equation 
cannot be argued by just the orthogonality. Since eigenfunctions evolve independently, 
the solution of (14.21) is 



M{r,x;q) 



2 — l—i^/2 roo 



K 2 q v/2 J n 



-u 2 t/S 



r 



V + %u 



Ki U (\/8qe x / 2 ) sinh(7ra) u du 



[-u/2] 



+ q u/2 1^ n \Y(-u-n + l) 

n=0 v ' 



( - 2n) e (-+ 2 -)V8 ir_^ 2n (v^e^ 2 ) . (4.4) 
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Therefore the option value can be written as 

„-Rt 



0(r,x) 



-vx/2-v 2 t/ 



2-nir j 



€+tOO 



dq e^ q M(T,x;q) . (4.5) 



Let us define the pricing kernel V(r, x, £) by 

»c+i oo 



V(t,x,0 



1 

2vri 



4 e (W2>-^ r/8 .M(r, x;q)dq 



<2 —^ I du sinh(Tra) « e -(« a +^Ws 



7T Z 



r 







x — / e^ q q- v/2 K iu (^qe xl2 )dq 



v + iu 



2tu 



-f/2] 



_j_ 2 3-^/2 e (l-v/2)x ^ 



-v - 2n) 



n=0 
2 /-e+ioo 



n! r(— z/ — n + 1) 



3 n(i/+n)r/2 



Then the option value becomes 



0(t,x) 



At e x 



V{r,x,Om^- 



(4.6) 



(4.7) 



(4i 



From this formula we can see that V(t, x, £) is essentially the probability density transition 
function of the stochastic process V T = J Q T e x dr. 



5 Integrating out q 



In this section we will work out the inverse Laplace transformation in (14. 7\i . For this we 
consider the following integration 



2tu 



/ e^q- u/2 K p {^~qe x/2 )dq. 

J t—ioo 



(5.1) 



To calculate this inverse Laplace transformation, we use the contour as in Figure [TJ It 
can be shown that the integration along tend to zero as A — > oo. When u+ |Re p\ < 2 
the integration along Cg also vanishes as 5 — > 0. Therefore the original integration along 
L is related to the integration along L + and L~ . Nevertheless it can be checked that 
the result we obtain in this way is still true for more general v. Note that due to the 
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Figure 1: The integration contour for the calculation of the inverse Laplace transformation 
(|5.ip . The radius of and Cs are A and 5, respectively. 



multivaluedness of the integrand, its values along L + and L are not same. We choose 
argg = 7r on L + , while argg = — ir on L~. Explicitly we have 



-as {S L . + L) eiqq "" 2K ' ( ^ qe " 2)iq 



1 

2vr7 



1 

2vr7 



(-dr) e-^ r {re^)~ u/2 K p {V^ r 1/2 e iw/2 ) 

POO 

+ / (-dr) e^ r (re- tn )- v/2 K p (V8^ ' r 1/2 e' i7v/2 ) 
Jo 

POO 

e -inu/2 / dre^ r r-^ 2 K p (V8^r l/2 e i7T/2 ) 
Jo 

POO 

-e 1 ™' 2 / dre-^r-^K^Vs^r^e^/ 2 ) 



1 ^("-l)/2 

2vr 



x I e- i7Tu/2 W^i p 

2 ' 2 



2-z/ + p^ ^ (2 - v - p 
2e x 



e in \ + e inu/2 



2e x 



(5.2) 
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In the last line above we have used the formula (6.643.3) of [20] to evaluate the integration 
in terms of the Whittaker function W K)ll (z). Note that z = is the branch point of 
W K . n(z), so W K)ll (ze m ) 7^ W KilJl (ze m ). To further simplify (15^2]) , we use the relation 
between W K ^(z) and the other Whittaker function M Kifi (z) 

W K Az) = F ^ ( " 2/X) , M Ki ,(z) + pn r(2 ^ , M K ^(z) , (5.3) 

r (± - k — fj,) r [± — K + fjL) 

and the Rummer's relation M K ^(z e ±ln ) = e ±tn (p-+ 1 / 2 ) M_ KitM (z) to take out the minus sign, 
together with T(l/2 + z) T(l/2 — z) = it/ cos(ttz), then we have 

/ = -tL= ^~ 1)/2 e~ e ^ Wi=, £ f ^ (5.4) 



2 ' 2 



Therefore the pricing kernel V(t;x,£) in ( 14. 71) is 

1 /"oo /q x\ (1— v)/2 



'2 e 1 

X Wl-v in 
2 ' 2 



r 



. \ 2 
Z/ + IU 



2 



n=0 



2(-„ - 2n) ( „ + „ )t/2 (2j?\ W f 2£ 

^ n!r(-i/-n + l) V f ' ~ 



sinh(7rw) u du (5.5) 



2 ' 2 



6 Integrating out £ 



In this section we will specify payoff functions for put and call options and then study 
their values. 

6.1 Put options 

The payoff function for the asian put option is 

<M0 = (Kt-Z) e{Kr-t), (6.1) 

where #(•) is the Heaviside step function. According to f)4.8p its value is 

-Rt r Kr 
At e x 



Op = I (KT-0V(T,x,$d£. (6.2) 
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By investigate the expression (15. 5p of the pricing kernel V(r, x, £), we see that we need to 
consider the following type integration 



Kt /o ,\(1-v)/2 







2e x V " ( e x \ „. /2e 



(lfr-0 ( T ) exp f -j 1 ^ I ^ ) di 



->X 



/oo 
(1 -y)y-W-W exp (-JL) ^ (J_) dy 

= 4e*(2k)We*p (-^) WL^ (£) , (6.3) 

where = KT/(4e x ), and we have use the formula (7.623.7) of [20] in the last line. 
Therefore the value of an Asian call option is 

Op = e —- e-^ + ^ T l\2k)^l 2 e- l l { ^W_^^(\ 
2tt 2 t L 2 ' 2 V 2k 



x 



r 



- ^ 2 

V + IU 



sinh(7ru) u du (6.4) 



gjrt- - 2 (-,/-2n) (, +n)r/2( )(3+ ,)/2 1/(4.)^ / 1 \ 

n=0 v ' v 7 

By using the relation (see (9.237.3) of [20] 0) 

W_3±^,_|_ n (2) = (-l) n n!^- n -( 3+ ^ 2 e- 2/2 L^- 2n (^) (6.5) 
with L~ u ~ 2k (z) being the generalized Laguerre polynomial, and 

W„-i/2,„(*) = z 1/2 -"e* /2 r(2p,z) (6.6) 

with T(2/i,z) the incomplete Gamma function, it can be shown that (16 .4p is exactly 
equal § to the result obtained in [17] through a different approach. [17J is based on an 
equivalence between two stochastic process, while our method seems more elementary, 
just doing integrations. Actually we can derive that equivalence by using the formulation 
in this paper. 

6.2 Call options 

Now we consider the Asian call option, whose payoff function is 

0c(O = U-Kt) OU-Kr). (6.7) 



2 The factor nl is missed in [20] . 

3 Need the replacement T here = 4T t hcrc and e x horo = So there- 
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The integration we need to do is 



AY 



2e x \ {1 ' u)/2 ( e x \ /2e a 



I c = / (£ - Kt) — exp -- Wi=», [-r)dZ 



= 4e 2 W 3+ ' )/2 /V-^- 3+ ^ rfy 

= 4e-(2^)(^exp (-^) WL^, f (£) • (6.8) 

We have used the formula (6.623.8) of [20] in the last line. Note that, unlike (16 .3p is always 
true, (I6.8p is convergent only under the condition v + |Rep| < —2. For the integration 
part of the pricing kernel V(r, x,£) to be convergent, we should require v < —2. Then 
there are at least n — 0,1 two terms in the finite summation part of (15. 5p . But these two 
terms are both divergent since u + (—v — 2n) > —2 for n = 0, 1. Therefore for call options 
we cannot naively interchange the order of the integrations over u and £. In the theory 
of finance we have so-called put-call parity relation [16] 

O c = Op + 1 ~ 6 RT e x - e~ RT K . (6.9) 
Kt 

We can use this relation to obtain the value of Asian call options from that of put options. 
However it is interesting to find a direct calculation method of this result. 

A The completeness proof 

In this appendix we will show that, for q > 0, 
1 f°° 

I(x,x') ■= — K iu (^e x/2 )K iu (^/8qe x ' /2 ) smh(7ru) udu = S(x - x') . (A.l) 
71 Jo 

Since the integrand is an even function due to K~i U = K^i U , we can extend the integration 
range to (—00,00). Therefore 

• riA 

I(x,x') = — lim / K„(J^e x/2 )K„(J8q^e x ' /2 )sm(iiv)vdv. (A.2) 

2tH A^ooJ_ iA 

By using the relation K v {£) = 2~ 1 n(I- u (z) — I v (z))/ sin(7rz/), we decompose I(x,x r ) into 
three terms 
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where z = y/8qe x '^ 2 and z' = y/8qe x '^ 2 . The first two terms are actually equal by 
interchanging v and —v. Since the integrands are holomorphic in the complex z/-plane, 
we can deform the integration path to a semicircle Ca := {Ae J<?i | — | < < |}. Therefore 



A^oo | sm(7rz/) 



" (/?/ + I_„(z) W)] \ . (A.4) 



When the order v is large and z is fixed, we have 



iv{z) ~ r(i + y) v2/ ' l ar s z l <7r - ( A - 5 ) 

By carefully analyzing the asymptotic behavior along C\ when A — > oo, it can be shown 
that |22j the first term of (1A.4j) actually tends to zero. So we have 



1 ' j 8 a^oo y CA sin(7rz/) I r(i + 1/) r(i -v) r(i - i/) r(i + v) 



2 



A^oo 8lT ir \ \z'J \z' 



A 00 f ( x—x' \n \ n I x—x' \n A n 

hm - y ( ~ + (-ir^T 

a^oo 87r ^—^ n! n\ 



n=0 



2 



In the second line we have used T(l + v) T(l — v) = nv/ sin (717/), and in the last line the 
limit representation of the Dirac delta function 

sin(Aw) . 
hm — i — - = 5(u) . (A.7) 

A^oo 7YU 

Therefore we have proved the completeness relation I 1A.1I) . 
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